21cm Absorption by Compact Hydrogen Disks Around 
Black Holes in Radio-Loud Nuclei of Galaxies 
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The clumpy maser disks observed in some galactic nuclei mark the outskirts of the accretion disk 
that fuels the central black hole and provide a potential site of nuclear star formation. Unfortunately, 
most of the gas in maser disks is currently not being probed; large maser gains favor paths that 
are characterized by a small velocity gradient and require rare edge-on orientations of the disk. 
Here we propose a method for mapping the atomic hydrogen distribution in nuclear disks through 
its 21cm absorption against the radio continuum glow around the central black hole. In NGC 
4258, the 21cm optical depth may approach unity for high angular-resolution (VLBI) imaging of 
coherent clumps which are dominated by thermal broadening and have the column density inferred 
from X-ray absorption data, ~ 10 23 cm" 2 . Spreading the 21cm absorption over the full rotation 
velocity width of the material in front of the narrow radio jets gives a mean optical depth of ~ 0.1. 
Spectroscopic searches for the 21cm absorption feature in other galaxies can be used to identify 
the large population of inclined gaseous disks which are not masing in our direction. Follow-up 
imaging of 21cm silhouettes of accelerating clumps within these disks can in turn be used to measure 
cosmological distances. 

PACS numbers: 98.54. Aj, 98.62. Js, 95.30. Jx, 98.80.-k 
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I. INTRODUCTION 

The discovery of an edge-on H2O maser disk in the 
nucleus of NGC 4258 [l], Q provides the best current ev- 
idence for a nuclear black hole outside the Milky- Way 
galaxy. The radial extent of the maser ring, ~ 0.14- 
0.28pc, and its measured Keplerian rotation profile be- 
tween ~ 770-1100 km s _1 , imply a central black hole 
mass Mk4x 10 7 M Q [3|. To obtain large maser gains, 
any such disk needs to be viewed at a nearly edge-on 
orientation with the maser clumps lying along the lines 
of minimal variation in the line-of-sight velocity 0, 0]- 
Given the small geometric likelihood of edge-on orien- 
tations, the subsequent discovery of maser disks in the 
nuclei of other galaxies d, 0, Q , demonstrated that cold 
nuclear disks must be ubiquitous. The maser clumps po- 
tentially mark the outskirts of the accretion disk that 
fuels the central black hole in these galaxies. 

The Milky Way galaxy does not show evidence for a 
compact nuclear disk but instead features young mas- 
sive stars at an orbital radius of < 0.3 pc from the cen- 
tral black hole, SgrA* [nj. These stars could not 
have formed in a central molecular cloud because of the 
strong tidal field of SgrA* [llj]. Instead, they are conjec- 
tured to have formed in a cold compact disk, at a sub- 
parsec distance from SgrA*, where the disk self-gravity 
was important T3, 13, 14|. For the generic radial profile 
of a viscous accretion disk, the Toomre Q-parameter is 
indeed expected to ap pro ach unity at the required dis- 
tance scale llj ill 17, 18 1, making the disk unstable for 
fragmentation. This coincidence could also explain the 
clumpy nature of maser disks at a rotation velocity of 
~ 10 3 km s _1 [lj|. Elsewhere in our Galaxy, masers are 
observed to be associated with regions in which massive 



stars form hJ 2(| . It is therefore likely that the observed 
maser disks in external galactic nuclei represent sites of 
star formation 14]. The abundance of maser disks im- 



plies that galactic nuclei experience common episodes of 
disk formation, during which massive stars - like those 
observed near SgrA*, are born. Winds from these mas- 
sive stars or supernova explosions (supplemented by ad- 
ditional feedback from the accreting nuclear black hole) 
could have dispersed the central disk in the Milky- Way 
galaxy on a short timescale of < 10 7 years. In the future, 
a new disk might form when fresh cold gas will assemble 
once again around SgrA*. 

Little is known observationally about the global gas 
distribution in the observed maser disks Here we 
propose to map atomic hydrogen in nuclear disks through 
its 21cm absorption of background radio glow around the 
black hole. It is known that radio-emitting plumes (which 
define the base of large-scale jets [2lj) exist around NGC 
4258 (see Fig. 1 in Ref. [H), but the angular extent of 
the coronal radio emission at 1.42GHz and its potential 
absorption by the disk have not been explored as of yet 
in the published literature In §2 we calculate the 

expected optical depth for 21cm absorption in nuclear 
gaseous disks. This absorption signature could be identi- 
fied through Very Long Baseline Interferometry (VLBI) 
observations of regions where the continuum backlighting 
is sufficiently bright. 

In principle, the 21cm absorption by nuclear disks can 
be mapped at high angular and spectral resolutions. The 
velocity and acceleration of clumps within the disk can 
then be used to infer the angular diameter distance to 
the sources, as demonstrated for maser clumps in NGC 
4258 @, HE |25|. The 21cm absorption feature can 



also be searched for spectroscopically (without spatial 
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resolution) in a survey over a large number of compact 
radio sources. Because the 21cm absorption signature 
would appear for arbitrary disk inclination, a dedicated 
search for this signature would be more likely to find 
nuclear disks than searches for masers which are limited 
to edge-on orientations of the disks. Follow-up VLBI 
imaging of disks could then be used to infer the central 
black hole mass for a large sample of galaxies. 



II. OPTICAL DEPTH OF COMPACT NUCLEAR 
DISKS 

The radiative transfer equation for the intensity I v of 
the 21cm line along a particular line-of-sight reads [2^ |. 



dl„ 
ds 
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[n 2 A 21 - (niB 12 - n 2 B 21 ) I v ] , (1) 



where v is the photon frequency ds is the path element, 
4>(y) is the line profile function normalized by J <j){v)dv = 
1 (with an amplitude of order the inverse of the frequency 
width of the line, Av), subscripts 1 and 2 denote the lower 
and upper levels of the line, n denotes the number density 
of atoms at the different levels, and A and B are the Ein- 
stein coefficients for the transition between these levels. 
We make use of the standard relations: B 21 — {gi/ g 2 )B 12 
and Bi 2 = {g 2 / gi)A 2 in/ I y , where g is the spin degener- 
acy factor of each state. For the 21cm transition, A 2 i — 
2.85 x lCT 11 ^- 1 and g 2 jg x = 3 [27]. The relative popu- 
lations of hydrogen atoms in the two spin states defines 
the so-called spin temperature, T s , through the relation, 
(712/711) = (g 2 /gx)exp{-E/kT s }, where E/k = 0.068K 
is the transition energy. In the regime of interest here, 
T s > E and so [(g 2 /gi)(ni/n 2 ) - 1] » E/kT s , and 
n 2 w jriff, where uh is the total number density of hy- 
drogen atoms. Moreover, the brightness of the sponta- 
neous 21cm emission is too weak to be detectable. We 
therefore ignore the first term in the square brackets of 
Eq. ((T|) and consider the absorption signature of the gas 
against a bright radio continuum glow in the background. 
Defining the optical depth along a ray as r = — Aln/„, 
we get 



. 3 h 3 c 2 A 21 N H 

T(iy) = 32^-^ [ ^ )] fclV 



(2) 



where Nh = / rinds is the column density of hydrogen. 

To get an estimate for the average optical depth value 
across the line profile, we write v(f)(v) = (Ai//i/) _1 , where 
(Av/vq) = (Au/c), is the fractional Doppler width of the 
line, corresponding to a velocity spread Av among the 
absorbing atoms. The minimum line width attainable is 
dictated by the spread in the thermal velocities of the 
atoms, for which [26|, Av = t>th = (2fcT/mjy) 1 / 2 , where 
m,H is the hydrogen atom mass. A larger width can be 
induced by a gradient in the bulk velocity of the gas 



along the line-of-sight and is calculated in the Sobolev 
approximation 37]. When the 21cm line is optically-thin 
(r < 1), the absorption signal obtains a width that re- 
flects all the contributing gas elements within the angular 
resolution and frequency band of the observations. 
Substituting all the coefficients into Eq. |2]) yields, 
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N H 



10 23 



x 10 3 K 



Av 



10 km s- 



(3) 

The specific numerical values that were substituted on 
the right-hand-side of Eq. ([3]), correspond to the ex- 
pected parameters of the atomic hydrogen disk in NGC 
4258 0] . X-ray observations of this system indicate a hy- 
drogen column density of ~ 10 23 cm -2 [3(J. Based on the 
observed X-ray luminosity of NGC 4258 and the warped 
geometry of its disk, the gas is expected to be predom- 
inantly atomic (rather than molecular) outside a radius 
~ 0.3 pc (see Fig. 23 in Ref. Q). Theoretical calcula- 
tions [H suggest that the atomic hydrogen [38| in the disk 
of NGC 4258 has an asymptotic temperature ~ 8000 K, 
providing a thermal velocity width of Vth ~ 10 km s _1 . 
At the high densities under consideration, we assume that 
the spin temperature T s is in collisional equilibrium with 
the kinetic temperature of the gas. High resolution imag- 
ing of nuclear disks can therefore be used to constrain 
their density and temperature distributions through Eq. 

The scale height of a thin accretion disk, h, is ex- 
pected [ll| to be a fraction ~ (vth/v<f>) of its radius r, 
where u^(r) = (GM/r) 1 ' 2 is the rotation speed at r for 
a black hole mass M. A line of sight which crosses the 
disk at an angle 9 relative to the normal to the disk sam- 
ples a spread in the line-of-sight bulk velocity that is 
< [2h/ cos9)[{dv <j) /dr)sm9} = (sin 2 9/cos9)v th . For the 
warped (bowl-shaped) maser region of NGC 4258, the 
value of cos# « 0.3 [3j yields Av ~ (1-3) x u t h- 

In order for VLBI imaging to reach the thermal broad- 
ening minimum of Av, a particular spatial resolution el- 
ement needs to be dominated by a single clump of gas 
with a coherent bulk velocity. Such a clump would have 
a small line-of-sight bulk velocity if it is located in front 
of the black hole and up to the full rotation speed on the 
side. If the clump fills only a fraction / of the source area 
within the resolution element, then r will be reduced by a 
factor of /. For the diffuse gas in the disk, the spatial res- 
olution required to achieve the thermal width minimum 
is of order the disk scale height, h ~ (vth/v^r. This 
resolution scale corresponds to ~ 10~ 2 r ~ 3 x 10~ 3 pc 
for NGC 4258. At a wavelength of 21cm and a source 
distance of 7.2Mpc, this resolution requires an unrealis- 
tic baseline of ~ 5 x 10 5 km, larger by a factor of 40 than 
the diameter of the Earth. Thus, a terrestrial VLBI will 
resolve the disk in NGC 4258 only around and outside 
the maser region (r > 0.14 pc). Coincidentally, this is 
indeed the region expected to be dominated by atomic 
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hydrogen An analogous disk around a quasar black 
hole whose mass M is larger by two orders of magnitude 
than in NGC 4258, could in principle be resolved out to 
a distance of ~ 1 Gpc. 

We conclude that the high value of the optical depth 
in Eq. applies to silhouettes of coherent clumps in 
which thermal broadening dictates the velocity spread 
Aw. Such clumps are expected to exist outside the maser 
region, where the Toomre-Q parameter is of order unity 
or lower 14|. If individual clumps of atomic hydrogen 



are not resolved or if the disk is smooth, then the optical 
depth would be diluted over a broader velocity width. In 
general, the absorption depth at a given frequency bin 
scales as the fractional (brightness-weighted) area of the 
continuum source over which hydrogen atoms resonate 
with photons in the observed frequency bin. 

Under a uniform background illumination, the absorp- 
tion spectrum of an unresolved circular disk can be ob- 
tained through a sum over concentric rings in the disk 
plane. We assume that the normal to the disk plane 
is inclined at an angle 9 relative to the line-of-sight. A 
single optically-thin ring with a circular rotation veloc- 
ity v<j,(r) — (GM /r) 1 / 2 and a radius r, gives a U-shaped 
spectral profile in terms of — 1 < 5(v) < 1 [32j |. 



(4) 



where S = [(u — Uq)/i/o]/[0.5Av/c], Av = 2v<p(r) sm9, and 
,- 2 ( N H (r) \ / T, 



T ri „g(r) = 0.44 x 10" 



V 10 23 cm- 2 / \8 x 10 3 K 



v^r) sin (9 
800 km s- 1 



(5) 



The total absorption feature of an unresolved disk can be 
obtained by summing over all the rings in which atomic 
hydrogen resides, weighted by the brightness distribution 
of the backlighting at vq = 1.42GHz. 

For an arbitrary background illumination, the net 
deficit in the fractional spectral intensity across the area 
S of an unresolved optically-thin source is given by, 



AI S 
Is 



f s I„(x,y)T(v,x,y)dxdy 
f s I„(x,y)dxdy 



(6) 



where (x, y) are the sky coordinates and the unabsorbed 
continuum source can be assumed to have a smooth (typ- 
ically power-law, I u oc v a ) spectrum across the absorp- 
tion line profile in the numerator. For a thin disk which 
is not perfectly edge-on, the exact expression for the op- 
tical depth t(u, x, y) in Eq. ([2]) can be approximated as 
the thermally broadened value in Eq. ([3]) at the Dopplcr 
shifted frequency ^o[l— t?ii/c], where v\\ (x, y) is the line-of- 
sight component of the bulk velocity of the gas. Clearly, 
in order for the spectral deficit to be noticeable, a dom- 
inant component of the radio emission needs to origi- 



nate behind the absorbing disk. If the background illu- 
mination originates from a narrowly collimated jet (as 
indicated by the 22GHz image of NGC 4258), then the 
absorption feature will be characterized by the low line- 
of-sight velocity spread (Av) of the material in front of 
the jet. In this case, the spectral deficit will be larger 
than the deficit associated with the full velocity spread 
of the disk. For the narrow jets of NGC 4258 we esti- 
mate Av ~ 100 km s _1 and r ~ 0.1. A dedicated search 
for the 21cm spectroscopic feature in other compact ra- 
dio sources can be used to identify new nuclear disks in 
distant galaxies. 



III. DISCUSSION 

The parameters of the maser disk in NGC 4258 imply 
that atomic hydrogen within the compact gaseous disks 
in galactic nuclei can produce measurable 21cm absorp- 
tion against the backlight of continuum radio emission 
around the central black hole. For NGC 4258, the 21cm 
optical depth may approach unity in silhouettes of co- 
herent clumps which are dominated by thermal broad- 
ening and have the column density inferred from X-ray 
absorption data 30] . Spreading the absorption across the 
rotation velocity width of the material in front of the col- 
limated jets in NGC 4258 results in an expected optical 
depth of - 0.1. 

More than half of all nuclear H2O megamasers show 
X-ray absorption with column densities Nh ~ 10 24 - 
10 25 cm" 2 (see Ref. ^ and Fig. 7 in Ref. Q), at 
which the optical depth for 21cm absorption might ex- 
ceed unity. High-resolution VLBI images of 21cm absorp- 
tion can be used to map the distributions of the density, 
temperature, and line-of-sight velocity of atomic hydro- 
gen in nuclear disks. Such maps could show direct evi- 
dence for spiral arms [3^], which are conjectured to exist 
based on the latest maser data in NGC 42580, EH . More 
generally, the maps hold the potential for testing current 
models of accretion disks, shedding light on the geometry 
of obscured (Compton-thick) quasars, and improving our 
understanding of star formation in galactic nuclei. 

A comprehensive search for a spectroscopic absorption 
feature in radio-loud nuclei of galaxies can be used to 
find a large number of inclined gaseous disks which are 
not masing in our direction. The improved statistics of 
known nuclear disks would provide better understanding 
of the duty cycle of black hole fueling and star formation 
in galactic nuclei. 

Remote atomic hydrogen within the host galaxy would 
also result in absorption but will be limited to low- 
velocity widths. The compact nuclear disk is expected to 
dominate the wings of the 21cm absorption profile which 
extend out to velocity offsets of ±10 3 km s _1 . Follow-up 
imaging of the nuclear disk can be used to separate out 
extended galactic absorption. 
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VLBI measurements of the velocity and acceleration 
of coherent hydrogen clumps within the disk can be used 
to infer the angular diameter distance to the source, as 
demonstrated with maser clumps @, [3, 24, [25|. De- 
tection of suitable radio sources out to sizeable redshifts 
could potentially place new constraints on the equation of 
state of the dark energy through the dependence of the 
angular diameter distance on source redshift 23|. The 
selection of a large number of suitable targets would be 
a particularly attractive goal [3ti| for the planned Square 
Kilometer Array [39i | . 



Acknowledgments. I thank Jim Moran and George 
Rybicki for useful discussions. 



[9 

[io; 
in 

[12 

[13 
[14 

[15] 



M. Miyoshi, J. Moran, J. Herrnstein, L. Greenhill, N. 
Nakai, P. Diamond, & M. Inoue, Nature 373, 127 (1995). 
J. Moran, L. Humphreys, L. Greenhill, M. Reid, & A. 
Argon, IAU Symposium 242, 391 (2007); ArXiv e-prints 
707, arXiv:0707. 10321 

J. R. Herrnstein, J. M. Moran, L. J. Greenhill, & A. S. 

Trotter, Astrophys. J. 629, 719 (2005). 

D. A. Neufeld, & P. R. Maloney, Astrophys. J. 447, L17 

(1995). 

D. A. Neufeld, P. R. Maloney, & S. Conger, Astrophys. 
J. 436, L127 (1994). 

L. J. Greenhill, et al, Astrophys. J. 582, Lll (2003); P. 
T. Kondratko, et al., Astrophys. J. 638, 100 (2006); P. T. 
Kondratko, L. J. Greenhill, & J. M. Moran, Astrophys. J. 
652 136 (2006). 

J. Braatz, et al. IAU Symposium 242, 402 (2007); 

J. Braatz, L. Greenhill, M. Reid, J. Condon, C. Henkel, 

& K. Y. Lo, IAU Symposium 242, 399 (2007); J. Braatz 

et al., IAU Symposium 242, 402 (2007) 

A. M. Ghez, in Massive Stars in Interactive Binaries, ASP 

Conf. Ser. 367, 689 (2007). 

R. Genzel, & V. Kar as, IAU Symposiu m 238, 173 (2007); 

ArXiv e-prints 704, larXiv:0704.T28T1 

M. Morris, Astrophys. J. 408, 496 (1993). 

Y. Levin, & A. M. Beloborodov, Astrophys. J. 590, L33 

(2003) . 

R. Genzel, et al., Astrophys. J. 594, 812 (2003). 

M. Milosavljevic, & A. Loeb, Astrophys. J. 604, L45 

(2004) . 

J. Goodman, Mon. Not. R. Astr. Soc. 339, 937 (2003). 



[16] J. Goodman, & J. C. Tan, Astrophys. J. 608, 108 (2004). 
[17] P. I. Kolykhalov, & R. A. Sunyaev, Sov. Astron. Lett. 6, 
357 (1980). 

[18] I. Shlosman, & M. C. Begelman Astrophys. J. 341, 685 
(1989). 

[19] M. J. Reid, & J. M. Moran, Ann. Rev. Astr. Astrophys. 
19, 231 (1981). 

[20] M. Elitzur, Astronomical Masers, Kluwer Academic Pub- 
lishers (1992). 

[21] G. Cecil, A. S. Wilson, & C. De Pree, Astrophys. J. 440, 
181 (1995). 

[22] P. T. Kondratko, L. J. Greenhill, & J. M. Moran, in 

preparation (2008). 
[23] L.J. Greenhill, New Astr. Re v. 48, 1079 ( 2004); L.J. 

Greenhill, ArXiv e-prints 708, larXiv:0708.0423l 
[24] E. Maoz, et al. Nature 401, 351 (1999). 
[25] E. M. L. Humphreys, M. J. Reid, L. J. Greenhill, J. M. 

Moran, & A. L. Argon, Astrophys. J. 672, 800 (2008). 
[26] G. B. Rybicki, & A. P. Lightman, Radiative Processes in 

Astrophysics, John Wiley & Sons (1979), pp. 30-33, 288. 
[27] G. B. Field, Astrophys. J. 129, 536 (1959). 
[28] V. P. Grinin, Astrophysics 44, 402 (2001). 
[29] R. Barkana, & A. Loeb, Astrophys. J. 624, L65 (2005). 
[30] A. Fruscione, et al. Astropphys. J., 624, 103 (2005). 
[31] S. L. Shapiro, & S. A. Teukolsky, Black Holes, White 

Dwarfs, and Neutron Stars: The Physics of Compact Ob- 
jects, John Wiley & Sons (1983), p. 436 
[32] G. B. Rybicki, & D. G. Hummer, Astrophys. J. 274, 380 

(1983). 

[33] G. Madejski, C. Done, P. T. Zycki, & L. Greenhill, As- 
trophys. J. 636, 75 (2006). 

[34] J. S. Zhang, C. Henkel, M. Kadler, L. J. Greenhill, N. 
Nagar, A. S. Wilson, & J. A. Braatz, Astr. & Astrophys. 
450, 933 (2006). 

[35] E. Maoz, Astrophys. J. 455, L131 (1995); E. Maoz, & C. 
F. McKee, Astrophys. J. 494, 218 (1998). 

[36] N. Kanekar, & F. H. Briggs, New Astron. Rev. 48, 1259 
(2004). 

[37] In the Sobolev approximation [28|], the term v<1>(v)Nh in 
Eq. ^ is replaced by cuh / {dv\\/ ds) where dv^/ds is the 
line-of-sight gradient of the line-of-sight velocity of the 
gas. In the cosmological context of 21cm absorption by a 
uniform intergalactic medium which encounters Hubble 
expansion, this gradient is simply the Hubble parameter 
(see, e.g., Ref. [29]). 

[38] The molecular gas in the disk is cooler (~ 10 s K) and 
could potentially be probed through other absorption 
lines. 

[39] http://www.skatelescope.org/ 



